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Abstract 

In this paper we propose a refinement of the heat kernel/zeta function treatment of kink quantum 
fluctuations in scalar field theory, further analyzing the existence and implications of a zero energy 
fluctuation mode. Improved understanding of the interplay between zero modes and the kink heat 
kernel expansion delivers asymptotic estimations of one-loop kink mass shifts with remarkably higher 
precision than previously obtained by means of the standard Gilkey-DeWitt heat kernel expansion. 

PACS: ll.15.Kc; 11.27.+d; ll.lO.Gh 

1 Introduction 

In 1961 Skyrme [i] proposed the following bold idea about the bridge between Classical and Quantum 
Field Theory: "If plane wave solutions of the field equations in the classical domain emerge as light mesons 
in the quantum domain, solitary wave solutions make the quantum leap, becoming heavy baryons in a 
non-linear quantum field theory" . This image from hadronic physics prompted the task of studying the 
quantum nature and properties of topological and non-topological classical non-linear waves with far- 
reaching consequences, ranging from Cosmology and Gauge Theories to Condensed Matter Physics. In 
1962, see [ ], Skyrme, together with Perring, found that the sine-Gordon equation with its roots in the 
Geometry of Surfaces of constant (negative) curvature offered an ideal playground in which Skyrme's 
idea was materialized analytically. In the mid seventies, Dashen, Hasslacher and Neveu addressed the 
investigation of one-loop kink fluctuations by developing the /i-expansion of the (1 + l)-dimensional 
(j)^ and sine-Gordon scalar field models. In the two seminal papers [3, 4], without and with fermionic 
fields, the authors encoded their findings in the evaluation of kink mass quantum shifts up to one- 
loop order, as well as the computation of the semi-classical/WKB corrections to the Bohr-Sommerfeld 
spectra of the sine-Gordon breathers. In the review reports [5] and [6] later developments on the subject 
achieved by several research groups in the USA and the EU between 1998 and 2006 are reported and 
summarized. This renaissance in the analysis of topological soliton fluctuations addressed mainly, but 
not only, supersymmetric versions of the same type of phenomena and/or generalizations to models with 
several scalar fields, also touching on the very delicate point of the influence of boundary conditions and 
their interplay with supersymmetry plus the choice of renormalization criteria. 

Here we shall concentrate on an specific method to deal with kink fluctuations. The one-loop kink 
mass shift induced by the kink quantum fluctuations is due to three contributions: 1) The kink zero point 
energy collecting the energy of the kink ground state where all the fluctuation modes are unoccupied. 
2) The analogous vacuum zero point energy that must be subtracted. 3) The energy induced by the 
one-loop renormalization mass counter-term on the kink background (measured with respect to the same 
effect on the vacuum). All these three terms are divergent and it was decided in [7] -supersymmetric 
kinks- and [8] -bosonic kinks- to regularize these quantities by means of the spectral zeta function of 
the kink and vacuum second-order fluctuation operators. In the sine-Gordon and (j)'^ models the kink 
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Hessian operators are Schrodinger operators of the Poschl- Teller type whose spectral problem is solvable. 
Therefore, one can exactly compute the kink spectral zeta function. 

In other interesting models the spectral information about the differential operators governing kink 
fluctuations is usually grossly insufficient for gaining analytical knowledge of the spectral zeta functions. 
The idea in the references quoted above was to rely on the determination of the spectral zeta function via 
the Mellin transform of the heat trace of the same operator and then use the asymptotic expansion of this 
latter spectral function to obtain information about the kink fluctuation asymptotics. With different stim- 
uli from Mathematics [9, 10] and Physics [11], interesting developments concerning the high-temperature 
expansion of the kernel of the generalized heat equation associated with a differential operator of the 
Laplace or Dirac type took appeared in the mid sixties. In particular, Gilkey, following classical works on 
the heat kernel proof of the index theorem, unveiled the meromorphic structure of the generalized zeta 
function, showing that the residua at the poles of the spectral zeta function, determined from the Seeley 
coefficients of the heat kernel expansion, are certain topological invariants/characteristic classes. In turn 
DeWitt used the heat kernel expansion to deal with quantum fields on curved backgrounds. During the 
last twenty- five years the use of heat kernel/zeta function methods in Quantum Field Theory has become 
a very important area of research: we explicitly mention the reports and textbooks [12, 13, 14, 1-5] for the 
reader to become acquainted with this important subject. During the last decade we (and our colleagues) 
profited from the theoretical machinery available by applying these tools at our disposal to express the 
one-loop kink mass shift as a truncated series in the Seeley coefficients of the kink-Hessian-heat function, 
see [8, 16, 17, 18, 19], in several (1 + l)-dimensional field theoretical models of a single real scalar field . 

Our goal in this paper is to address a very delicate issue. The Gilkey-DeWitt heat kernel expansion 
is an "inverse temperature" power series such that the zeroth order term is the heat kernel of the "free " 
Laplace/Helmotz operator. Thus, the GDW heat kernel expansion coincides at infinite temperature with 
the free heat kernel. Nevertheless, because of the completeness of the eigenfunctions of the kink-Hessian 
Sturm-Liouville problem the corresponding heat kernel cannot go to zero at low temperature if there are 
zero modes in the spectrum, as happens in the GDW expansion. Given that there is always a zero energy 
kink fluctuation, the translational mode, here we propose to modify the GDW heat kernel expansion 
to incorporate this notion when heat trace/zeta functions techniques are to be applied to the study of 
kink fluctuations. Even though the standard lore is that zero modes enter at two-loop order in the 
/i-expansion of quantum mass shifts, they play a hidden role noticed, for instance, in the Cahill-Comtet- 
Glauber computation [20] in terms of the bound state eigen-values and the threshold of the continuous 
spectrum. Thus, we shall assign special status to the zero modes that, somehow, is tantamount to 
performing the heat kernel expansion around the sine-Gordon kink Hessian. The new procedure is not 
only conceptually more satisfactory but enhances the numerical precision in the computation of the kink 
mass quantum correction to a remarkable extent. 

The organization of this paper is as follows: in Section §.2 we briefly review the semi-classical kink mass 
DHN formula and variations of this expression obtained in the heat kernel/zeta function regularization 
procedure. Section §.3 contains the main novelty in this paper: the modified heat kernel expansion 
adapted to the existence of the kink zero mode is introduced. In the Section §.4, the new formula is 
applied to the sine-Gordon and (j)'^ kinks to test the modified procedure. The response found in these 
models is in complete agreement with the exact results, improving the approximation reached in the 
traditional heat kernel expansion approach. We shall also work this method to evaluate the one-loop 
mass shifts in the scalar field models recently addressed in [19]. The interest in these models lies in 
the fact that full information about the kink Hessian spectrum allows us to know exactly the quantum 
correction offering a very good playground to gauge all the methods. Finally, we shall run the new 
procedure in two one- parametric families of scalar field models recently analyzed in [ J n] : a generalized 
family of 0^ models and the double sine-Gordon model. In both cases the DHN formula is not applicable 
and the exact result is not known. Nevertheless, interesting results can be obtained in the heat kernel 
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approach, old in [18] and new in this paper. In the Appendix, a Mathematica code that automatizes the 
computations involved in the new procedure has been included. 



2 One-loop quantum kink fluctuations 

In this Section we briefly review the standard lore about the conceptual understanding of quantum 
kink fluctuations as well as the approach of our group to this subject in the heat kernel/zeta function 
framework. We shall follow the notation and conventions fixed in the recent reference [18]. 



2.1 Classical field theoretical models and kinks 



The action governing the dynamics in a (1 + l)-dimensional relativistic one-scalar field theoretical model 
is of the form: 



Here, tp{y^) : M ' — )■ M is a real scalar field; i.e., a continuous map from the (l-l-l)-dimensional Minkowski 
space-time to the field of the real numbers, ip = t and = y are local coordinates in M^'^, which is 
equipped with a metric tensor that we choose: gfj_y = diag(l, —1), /x, z/ = 0, 1. 

We shall work in a system of units where the speed of light is set to one, c = 1, but we shall keep 
the Planck constant h explicit because we shall search for one-loop corrections, proportional to h, to the 
classical kink masses. In this system, the physical dimensions of fields and parameters are: 



[h] = [S] = ML , [y^]=L 



M5L5 



[il] = ML-^ . 



The models that we shall consider are distinguished by different choices of the part of the potential 
energy density that is independent of the field spatial derivatives: f7[^(y'^)]. In all of them, there will 
be two special parameters, and 7^, to be determined in each case, carrying the physical dimensions: 
[m^] = and [7,^] = M~2L~2. We define the non-dimensional coordinates, fields and potential in 
terms of these parameters: 



mdy^i , xq = t , xi = X , (/) = 7dV 



The action and the "static " part of the energy are also proportional to dimensionless action and energy 
functionals, namely: 



^S[ct>] 



ill I 



Td 



dx^dx^ 
dx 



1 



2 dx^ dxf^ 



uicpixn] 



(1) 

(2) 



where we shall assume that U{(j)) is a non-negative twice-differentiable function of (p: U{cl)) G C^(M) 
and U{(j)) > for (p £ M. The configuration space C of the system is the set of field configurations 
C = {(t>{to,x) e Maps(M\M)/£;[(/)] < +00}. We assume that the set M = {(j)^'^ / U{(t)^^) = 0} of the 
minima of C/ is a discrete set. The simplest solutions of the field equations 



(f>{t,x) 



6U 



it,x) 



(3) 
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are static and homogeneous, and hence the elements of A4. The small (quadratic) fluctuations around 
any of these constant solutions, (j){t,x) = + 5(j){t,x), satisfy the linearized field equations: 



w^-^2+A^^^t^^)+<^mf] = ^ with 0" - -^ 



. (4) 



The solutions of (4), of the form 6(j)k{t,x) = e*'^*-'^^*/fc(x), are the normal modes of fluctuation of the 
system near one minimum of U. These linear waves are built from the eigenfunctions of the second-order 
"vacuum"^ fluctuation differential operator: 

Ko = -^+v^ , Kofkix) = u\k)fk{x) , fk{x) = e"'^ ,k€R , u\k) = + . (5) 

We remark on a point usually unnoticed: there is a "half-bound" state at the threshold of the continuous 
spectrum with constant eigenfunction. 

If the cardinal of 7W is greater than 1, there may exist spatially dependent static solutions of the 
solitary wave type (non-dispersive non- linear waves), which in the literature are referred to as kinks or 
lumps, see e.g. [21, 22, 23]. The static kink solutions 4>k{x) satisfy the first-order ODE ^ 



j^ = ±^j2im , (6) 

together with the asymptotic conditions guaranteeing finiteness of the energy: 

lim (l)K{t,x) = e 7W , lim 0i^(t,x) = (j)^^^^ £M , lim ^'^^j^'^^ = . (7) 

x— >+oo x— >— oo x— >iboo ox 

Small fluctuations around any of the equivalent kink/antikink solutions, (p{t,x) = 4>k{x) + (5(/)(t, x), still 
solving (3) satisfy the linearized field equations: 



^,-^,+v' + V{x)\6<t^{t,x) + Omf]=^ ' H 



= v'^ + V{x) . (8) 

</)x(a;) 



Again the solutions of (8), of the form 5(j)k{t, x) = e'^'^^'^^^ fq{x) , are the normal modes of kink fluctuations 
coming from the eigenfunctions of the second order kink fluctuation differential operator: 

dx"^ 



^ = -£2+^' + ^(^) • (9) 



We shall assume that the potential well goes to the same limit at both ends of the straight line: 
lima;-s.±oo ^(a:) = 0. This is true in all models for which = v"^ = where (/>^*^ and 

are the vacuum solutions that are connected by the kink. The kink normal modes built from the 
eigenfunctions of K linear superposition 

Kf,{x)=u\q)fg{x) , co^{q)=q^+v^ 

are not plane waves but some dispersive wave functions distorted by the kink. Contrary to Kq, which is 
a Laplace/Helmotz operator, X is a Schrodinger operator for a potential well. Therefore, the spectrum 
of K is formed by scattering states, bound states and, possibly, half-bound states. 

The lowest K-bound state fluctuation (mesons trapped by the kink) is always a zero mode: fo{x) = 
Wo = 0. This i^-eigenfunction is no more than the Goldstone boson due to the spontaneous symmetry 
breaking by the kink of the invariance of the system with respect to spatial translations: a; — t- x + a. 



^ We pass to use QFT terminology: the linear waves in the quantized theory are the light mesons, whereas the minima 
of U are the vacua of the system. 

^The kinks are also static solutions of the PDE equation (3). 
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2.2 Kink mass quantum correction induced by one-loop fluctuations 
2.2.1 The "first" DHN formula 

The problem of computing the shift in the kink mass induced by kink fluctuations in the one-loop order 
of the sine-Gordon and A(0^)i+i models was solved by Dashen, Hasslacher and Neveu in [•>]. They 
succeeded in building a formula that collects all the effects contributing to this quantum effect.^ We 
describe the formula as a regularization-renormalization construction in two steps. 



1. Zero-point kink energy renormalization: mode- by- mode subtraction of the zero point vacuum energy 

In the first step, a zero point kink energy (the energy due to all the -fC-fluctuation modes being 
unoccupied) renormalization is performed by means of a mode-by-mode subtraction of the zero 
point vacuum energy -all the Xo-fluctuation modes being unoccupied-. In the Appendix of [iT], we 
generalized the DHN formula to kinks other than the sine-Gordon and A(i?i>^)i+i kinks, using the 
one-dimensional Levinson theorem, finding the formula: 



lim r dk - ^Vk^ + v^ + :^{V{x)) + V + Sfecu, - ^1 , (10) 



2 La 

where {V{x)) = dxV{x) and d{k) are the phase shifts associated with the second order kink 
fluctuation operator (9). The two last terms in (10) collect the contribution of the possible half- 
bound state of K and the ever-present half-bound state of Kq and represent the novelty in [' ~] with 
respect to the DHN formula. 

2. Mass renormalization counter-term contribution: 

AEi{<j)K) in (10) is still divergent because in models with interactions there are more ultraviolet 
divergences than the vacuum energy. In (1 + l)-dimensional scalar field theory normal ordering 
tames all the uv divergences. The ordering induces a self-energy counter-term that, in turn, via the 
expectation value of the scalar field operator at the kink and vacuum coherent states, produces the 
following contribution to the one-loop kink mass shift expression, see e.g. [6]: 

AE2{cPk) = -^{V{x)) Ihn JL^ , (11) 

The addition of these two quantities, previously regularized separately by choosing a cutoff that counts 
the same number of eigen-modes over the kink as over the vacuum, provides the finite value of the 
one- loop kink mass shift: 

AE{(I)k) = AEi{(t>K) + AE2{^k) ■ 



2.2.2 The zeta function regularization procedure 

In order to apply the above DHN scheme effectively in the computation of one-loop kink masses, it 
is necessary to know the eigenvalues of the bound states and the scattering wave phase shifts of the 
operator K. Only for the sine-Gordon and kink is this information fully available and other 

approaches to the problem are necessary in the investigation of quantum fluctuations of more complex 
kinks. The application of the zeta function regularization procedure has proved to be very effective in 
such a situation. The vacuum energy induced by quantum fluctuations is flrst regularized by assigning 

^We refer to this formula as the first DHN formula, to be distinguished from the "second" DHN formula applicable to 
the quantum sine-Gordon breathers. 
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it the value of the spectral zeta function of the ii'o-operator (a meromorphic function) at a regular point 
in s E C: 

Ai?(0«) = ^Ci.o(-^) ^ Ai^(c/>»)H = (^XcKois) , (12) 

Z Z 771(1 K^d^ 

where /i is a parameter of dimensions L^^ introduced to keep the dimensions of the regularized energy 
right. We stress that s = — | is a pole of this function. The same rule is applied to control the kink 
energy divergences by means of the spectral zeta function of K and, finally, the kink Casimir energy is 
regularized in the form: 

AE,{c^k)[s] = AEo{^k)[s]-AEo{<P^'^)[s] = ^(^) ''\cj,{s) - Ck,{s)) . (13) 

The advantage of this procedure is that there is no need for detailed information about the spectrum 
of K to calculate the spectral zeta function. It is possible to elude this problem in two steps: 

1. First, an asymptotic formula for the i^-heat trace can be obtained from the GDW heat-kernel 
expansion in terms of the so-called Seeley coefficients. 

2. The Mellin transform of the heat trace in turn provides the K-zeta function such that the regularized 
kink energy reads: 

A E,{^k)[s] = ^ (^)'^' (^K{f3) - HkM)] ■ (14) 

The energy due to the one- loop mass counter-term (11) can be also regularized by the zeta function 
procedure: 

AE2{^kM = ^{V{x))(^X ' lim^^^^^Ci.„(. + l) . (15) 
Finally, we write the zeta function regularized DHN formula: 

AE{(I)k)= lim AEi{^k)[s]+ lim AE2{(I)k)[s] ■ (16) 



3 Zero modes and the heat kernel asymptotic expansion 

3.1 Conflict at low temperature between zero modes and the heat kernel factoriza- 
tion 

The spectral if-heat trace hxifi) = Tr^2 e~^^ admits an integral kernel representation in terms of the 
fundamental solution of the ii'-heat equation: 

hK{/3) = j^dxKK{x,x;(3) , (^-^ + KK{x,y; P) = , KK{x,y;0) = 5{x - y) . 

The trace, of course, is in the functional sense, $7 is a finite (but very large) normalization interval on 
the real line, and /? is a fictitious inverse temperature. From Sturm-Liouville theory, the eigen-function 
expansion of the heat kernel 

b 

KK{x,y;(3) = f^{y)fo{x) + J2rn{y)fn{x)e-^''" + dk fl{y) h{x) e-^^'^'^^ (17) 

n=l 
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is founded in solid grounds. Writing (17) we had in mind the general structure of the spectrum of the 
Schrodinger operator K governing the kink fluctuations: one zero mode fo{x), a discrete set of bound 
state eigenfunctions fn{x), and the continuous spectrum eigenfunctions fk{x) characterized by a "wave 
number" G M. 

In the low temperature /3 — t- +oo regime, only the contribution of the eigenfunction in the algebraic 
kernel of K survives, whereas the completeness of the eigenfunction expansion fixes the K-heat kernel at 
high-temperature /3 — >• 0: 

hm KK{x,y;(3) = fS{y)fo{x) , lim y; /3) = 5(rE - y) . (18) 

The standard GDW strategy (and many others) to find the solution of the parabolic equation 

KKix,y;l3) = (19) 



9 2 T./ ^ 



with a (5-source on the diagonal of M? at infinite temperature Kk{x, y; 0) = 6{x — y) is to start from the 
well known KQ-heaX, kernel 

KK,{x,y-p) = ^^e-^^\-'^ , fl{x) = ^e'^- , u:l{k) = + v"" , G M , 

and search for a iC-heat kernel in the factorized form: 

KK{x,y,P) = KK,{x,y;P)A{x,y-p) , A(x,y;0) = l . (20) 
The GDW heat kernel expansion 

Kk{x, y; /3) = -r^^'^'' ^~~Y1 y)/^" (21) 

^ n=0 

trades the PDE (19) by a recurrence relation between the Seeley densities an{x,y). Setting ao{x,y) = 1, 
one can prove, see e.g. [10], that the GDW heat-kernel expansion (21) is an asymptotic solution to the 
PDE equation (19), i.e., denoting the partial sums in the form K^^{x^ x; j3) = Kko{x, x; j3) Yln=o ^n{x, x)P^, 
we have: 



9 2 ^ 



K^'{x,x-p) = ^=R^^{x,x-^)P^^ , E:^*(x,x;0) = "5(0)" , (22) 



where the remainder R^^{x, x; /3) is a C°° function of x. 
Because 

lim KKo{x,y;P) = , lim KKo{x,y;P) = S{x -y) , (23) 

it is clear that although the asymptotic of Kx{x,y; P) read from (20) and (23) fits with (18) at high- 
temperature, it fails at low temperature due to the zero mode. This discrepancy forces us to limit the 
integration range in the Mellin transform to a finite one if we use the standard factorization. In other 
fields, e.g. heat kernel proofs of the index theorem, computation of anomalies in QFT, this is not a 
difficulty because these effects are related to residua at the poles of the spectral zeta function and the 
entire parts can be neglected. However, here we are in trouble. The main reason is that there is a value 
of /3 between low and high temperatures where the heat trace starts to depart from the exact one due 
to the zero mode. This causes many computational inaccuracies. A considerable amount of ingenuity is 
needed to obtain good estimations of the one-loop kink mass shift from the kink zeta function via the 
Mellin transform of the kink heat kernel expansion. 
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As a loophole, we propose the improved factorization 



Kk{x, y- P) = Kk, {x, y; (3) C{x, y; (3) + g{f3)e- — f^{y)fo{x) , 



(24) 



matching the asymptotic behavior of the ET-heat equation kernel at both low and high temperature, (18) 
and (23), provided that the /3 — t- oo and /3 — t- limit of C{x, y; j3 and g{l3) will be such that: 



lim C(x, -y; 0) = 1 



hm 5(/3) = 1 , hm g{P) = 
/3— s>oo /3— >0 



(25) 



Having modified the standard factorization to adapt the formalism to the existence of zero modes, we 
follow the GDW route. Plugging the ansatz (24) into the iT-heat equation (19) the following "transfer" 
equation for C{x,y;(3) arises: 



+^e^-'fSiy) 



dgiP)., ^^5(/3), . Mx) 



. 



(26) 



3.2 The modified asymptotic expansion of the kink heat function 

This PDE equation is traditionally solved by means of a power series expansion: 



C{x, y]l3) = ^ c„(x, y) , co(x, y) = I 



(27) 



n=0 



Note that co{x, y) = 1 is obligatory to comply with the infinite temperature value of C{x, y; /3). Plugging 
the power expansion in (26) this PDE is transformed into the recurrence relations: 



n=0 



, ,N . X d'^Cn{x,y) c„+i(x,y) , ^ . . 
(n + l)Cn+i{x,y) ^ h (x - y) — h V{x)Cn{x,y) 



dx 



/3" + 



+V47r/3e^'^Vo(y) 



. 



(28) 



We stress that the modified factorization ansatz leads to new recurrence relations as compared with 
those encountered by Gilkey and DeWitt: namely, all the terms in the second line of (28) are new. This 
structure suggests the optimum choice oi g{/3). Setting 



(29) 



we simplify maximally (28) in the sense that derivatives g{f3) do not enter. Moreover, the error function 
is selected as the best option to comply simultaneously with (29) and the asymptotic conditions (25). 

Prom the asymptotic expansion of the error function Erf z = -^e~^^ Y1^=q (2n+i)!! ■^^"'^"^i the expres- 
sion (28) leads to the recurrence relations: 



(n + 1) c„+i(x, y) - ^ +{x- y) ^^"^^^^'^^ + V{x)cn{x, y) + 



. 



(30) 
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The modified heat kernel expansion takes the form: 

KK{x,y;p) = KK,{x,y;P)Y,Cn{x,y)P'' + ^vme-^ r^{y)h{x) , (31) 

n=0 

where the Cn{x,y) densities satisfy the relations (30), whereas the diagonal kernel reads: 
Kk{x,x,I3) = hm KK{x,y, (3) = -j== > c„(x, x) /?"■ + Erf(/3)/o (x)/o(x) , c„(x,x) = limcn(x,y) 

^ n=0 

and the heat trace is 

/ix(/3) = / dxi^;^(x,a;,/3) = -== Vc„(K)/3" + Erf(i;v^) , Cn{K) = dx Cn{x , x) . 
Jn v47r/3 ^ Jq 

Thus, to determine the i^-heat function hxif^) is tantamount to identifying the Seeley coefficients Cn{K). 
This task is accomplished by solving recursively (30) but to do this it is necessary to deal with the following 
subtlety: the operations of taking the y — )• x limit and the derivatives with respect to x in the formula 
(30) do not commute. To cope with this problem we introduce the following notation: 

(^)C„(x) = hm ^^^^i^ , (32) 

y^x Ox'' 

the Seeley densities, for instance, b cing: Cfiix^x^ — ^^^C^(ir). MorGovcr, th.6 "initis-l" conditions fix tliG 
first coefficient ^^^C{}{x) , Vfc to be: 

i^)C,{x) = \i^pl = 5^^ . (33) 
y^x ijx'^ 

Taking the fc-th derivative of the recurrence relation (30) with respect to x and later passing to the y — >• x 
limit in the resulting recurrence relations we obtain: 



n + k 



('=-^)c„_,(x)-5:^^^)U(-^)c„_,(x) 



3=0 

-2.,/„(.)!|!p*„,„_,-A,.)^|^i^(I + 2.)] . (34) 

which affords us the identification of ^^^Cn{x) in a recursive way from the initial ones (33). Note that 
the computation of the density c„(x,x) = ^^^Cn{x) requires the coefficients (°)C„_i(x), (i)C„_i(x), 
(''+'^)Cn-i{x), which are in turn determined from the densities ^^^Cn-2{x), ^^^Cn-2{x), ■ ■ ■ ^'^'^^^ C n-2{x) , 
etcetera, until we reach the densities ^'^^Co(x) = 6^^. 

We list the first three densities c„(x,x) derived from (34) for the kink fluctuation operator K 

co(x,x) = (o)Co(x) = 1, 

ci(x,x) = WCi(x) = -V{x)-Avf^{x) , (35) 
C2(x,x) = (0)C2(x) = -^-^+^-(V{x)f + ^v^f,{x)+Avf^{x)V{x), 

as well as the first three Seeley coefficients 
co{K) = I, 

ci{K) = -{V{x))-^v, (36) 
C2{K) = -^-{v"{x)) + l{{V{x)f) + ^v' + 4v{V{x)fS{x)) , 
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where I is the length of the interval fl. 

Subtraction of the i^o-tieat trace from the K-heat trace expansion amounts to dropping the co{K) 
coefficient, and we find: 



hKi/3) = hK{f3) - hKoiP) 



^c„(i^)/3"-i+Erf(^;v^) . 



(37) 



n=l 



In practical calculations we shall truncate the hxiP) series to a finite number of terms^, which defines 
the function 



^c„(ir)r-5+Erf(^;v^) . 



(35 



n=l 



3.2.1 An example: the asymptotics of the Xcj)'^ kink 

To illustrate the power of this modified procedure in the calculation of one-loop kink masses, let us 
focus on the X{(j))l_^_l model. K is the Schrodinger operator for the second member of the hierarchy of 
"transparent" Poschl- Teller potentials, see [ . ]. The spectral problem is completely solvable and the 
-fC-heat trace, a function of /3, is known exactly : 



hKiP) = %(/3) - HkM = Erf ( + Erf (2^/3) , 



(39) 



see [8, 18]. In Figure la, the exact hxiP) "renormalized" heat function has been plotted. A comparison 
with the approximated heat function given by the truncated asymptotic series hK{/3,Nt) is shown in 
Figure lb for several values of Aj. We emphasize that we used the modified GDW heat kernel expansion. 
The Seeley coefficients needed to compute hK{f3,Nt) are listed in Table 2 of Section §.4. The agreement 
with the exact result is remarkable even for low orders of At. In Figure Ic different plots of hxiP', Nt), 




Figure 1: The exact function /i_ft'(/3) (left), partial sums hK{l3,Nt) from the modified factorization method 



(middle) and partial sums hxil^-, Nt) from the standard GDW factorization method (right) for the X{4>)'\ 



kink. 

the truncated heat functions derived from the GDW factorization (20), are shown as colored lines for 
increasing values of A^ from the bottom (red) to the top (violet): 



hK{P;Nt 



2 Nt 



(40) 



n=l 



The Seeley coefficients an{K) for the Xcj)'^ kink are taken from the References [8, 18]. The main point is 
perfectly clear here: 



lim hK^ = 1 = lim hK{P,Nt) , VA* 

/3— >-oo l3^oD 



lim hKW,Nt) = 



After all, it is an asymptotic series such that there is an optimum truncation order approximation to the exact value. 
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The exact heat trace and the asymptotic approximation through the modified GDW factorization for the 
A(/)'^ kink tend to unity when the temperature decreases to zero: /3 = oo. The asymptotic approximation 
from the standard GDW factorization, however, tends to zero at low temperature. It is clear that the zero 
mode is not captured in this latter case and one must choose an optimum value of /3 as the upper limit 
of the integral in the Mellin transform. There is a second point to be emphasized: hxil^jNt) approaches 
hxif^) much faster than hxiP, Nt). Higher Nt is necessary to fit in with the exact function in this 
latter case. Because Ck{s) — Cko{s) = ffsj Io° dfi l3^~^hK{P), it is evident that a better approximation 
to the "renormalized" kink zeta function is C,k{s) — Qko{s) = d(3 (3^~^hK{f3; Nt) rather than 

Ck{s) - Cko{s) = YTs) /o" d/3f3'-'hK{f^,Nt). 



3.3 Kink mass quantum corrections in one-component scalar field theory 

The Mellin transform of the renormalized K-heat trace expansion hxiP) (37) gives us the renormalized 
K-zeta function as the series (involving Euler Gamma functions) 



Cxis) - Cko{s) 



n=l ^ ' * 



sT[s\ 



which in turn, recall (13), prompts the following expression for the regularized kink Casimir energy 
^Ei{(I)k)[s]: 



^Ei{(t>K)[s] 



1 (y(x))r[. + ^] 
r[s] 



_2 r[g + ^] 



+ 



+ 



1 ^ Cn{K) V\s + n-\\ 

2-^ ,,2n+2s-l 
ra=2 



r.s 



Here we have used the explicit expression of c\{K) given in (36) to write the contribution of the first 
term of the series. Because the contribution of the regularized one-loop mass counter-term A£^2(i?^_ft')[s], 
see (15), is 



(n^))r[g + ^] 1 

2s+l 



a crucial cancelation, obeying the heat kernel renormalization criterion, see e.g. [7], occurs after the 
addition of these two regularized quantities and we obtain the following series for the regularized one- 
loop kink mass shift A^(0x)H = I^Ex{(\>k){s\ + AS2(0x)[s]: 



^E{^kM \ ( II 



2 \ml 



1 



'A-K 



n=2 



Cn{K) V[s + n-\] 



T\s + \] 



2n+2s-l 



T[s\ 



TT V 



2s 



1 r[g + i] 

^y2s sT[s] 



Still, the last two terms are divergent at the physical point s = which is a pole of T{s + 5). There 
is, however, a miraculous cancelation that is observed by taking the dangerous limit carefully: 



lim 
V (\ 



1 r[s + ii 



■n \e 



-[7 + 21ogi; + ?/;(-i)] + 2) + 02{e 



lim 



-V fl 
-TT \£ 

2v 

TT 



[j + 2logv + i;{-l)]] +Oi{e) 
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Thus, zero point and mass renormahzation get rid of all the divergences and we end with the renormalized 
one-loop mass shift formula derived from the modified asymptotic series expansion of the K-heat function: 

^^^-'-^t-^mi^'y-nn-l]. (41) 

We have elaborated a Mathematica code, shown in the Appendix of this paper, which automatizes the 
computation of the quantum correction to the kink mass in one-component scalar field theory models. 
The algorithm is based on the evaluation of the Nf-th order partial sum of the series (41), 

^^|f^ = -;-^X:^»W("V-"r|n-ll (42) 



n= 



as a good estimation of the quantum correction. This Mathematica program uses the recurrence relations 
(34) and the definition of the Seeley coefficients to figure out the coefficients Cn{K) in (42). The inputs 
needed in this program only include the potential term U {(p), which characterizes in action (1) the specific 
model chosen, the two vacua which are connected by the kink solution whose mass quantum correction 
we are interested in, and the order of truncation Nf to be considered in our estimations (42). 



4 Specific (1 + 1)-D one-component scalar field theory models 

In this Section we shall apply formula (42), in the environment of the Mathematica code, to compute the 
one-loop kink mass shift in several well selected (l-l-l)-dimensional scalar field theory models. We start 
with the very well known sine-Gordon and t/)^ kinks as a good test of the new procedure. We shall see 
that in these paradigmatic cases the new method works extremely well whereas the analysis helps in the 
conceptual understanding of the modified heat kernel expansion. We shall then move on to explore more 
exotic avenues: First, we shall apply the procedure to compute the kink mass quantum correction to the 
a = N "parent" models, see [10], a very favorable arena because the second-order fluctuation operator 
K = N"^ — -A^ S N belongs to the transparent Posch- Teller hierarchy. Second, a one-parametric 

family of generalized (f)^ models having kinks and "double" kinks, see [IN], will be also addressed and 
the one-loop kink mass shifts estimated. Here, to calibrate the quality of the results only analyze of 
self-consistency are available. The K operator does not belong to any class of Schodinger operators 
with known spectra. Therefore, no other procedures of computation such as the DHN or the Comtet- 
Cahill-Glauber formula are applicable. Third, the kinks in the double sine-Gordon model [ ] will be also 
considered. Again, we shall obtain estimations of the kink mass quantum shifts that are not accessible 
with other procedures. 



4.1 The sine-Gordon model 

The part of the energy density potential independent of the field spatial derivatives is in the sine-Gordon 
model: 

U{(t)) = 1 - cos (/) . 

The sine-Gordon kink (pKix) = 4arctane^' -|- 27rn, n E Z, is a soliton solution that connects the vacua 
27rn and 27r(n -|- 1). The second-order small vacuum and kink fiuctuation operators are respectively: 

Kq = -—^ + 1 , K = -—^ + l-2sech^x . 

Thus, V = 1 and V{x) = — 2sech^x, in agreement with the notation used in the previous Sections. The 
spatial derivative of the kink solution is the (normalized) zero mode of fiuctuation: /o(x) = 4=sechx. 



12 



From (34) and spatial integration of the Seeley densities we find that all the Seeley coefficients for the 
sine-Gordon kink vanish. This reveals the hidden nature of the modified heat kernel expansion: the 
choice g{f3) = Erfy^ means that the new expansion is with respect to the heat kernel of the sine-Gordon 
kink, not of a constant background!!. Therefore, = Erf(-^/^) is the exact sine-Gordon kink heat 

trace. The numerical estimation of the sine-Gordon kink mass quantum correction through (42) thus 
provides the exact result: ^^S'^'^^ = — ^ = —0.31831. The optimum number of terms in the asymptotic 
formula is A^f = !! in the always surprising sine-Gordon model. 



4.2 The A(0)^+i model 

The potential term in this model is: 



The kink solitary wave 4>{x) = tanhx connects the two vacua 0'-^^ 
kink Hessian operators are respectively: 



-1 and (/)(^) = 1. The vacuum and 



+ 4 



K 



+ 4 — 6 sech X 



such that V = 2 and V{x) = — Gsech'^rc. The zero mode is the spatial derivative of the kink, i.e., properly 
normalized, fo{x) = ^ sech^ x. The symbolic resolution of the recurrence relations (34) for the analytic 
expressions of V{x) and fo{x) in the Mathematica code provides Seeley densities that, integrated over the 
whole real line, become the Seeley coefficients entered into the formula (42). These coefficients are listed 
in the Table 1. In this table we have also specified the values of the one- loop mass shift obtained from the 
modified asymptotic series (42) for different values of the truncation order, A't, graphically represented 
by the blue line of the attached figure. For the sake of the comparison we also collect from [8, 18] the 
same data obtained from the standard asymptotic series, represented by the red line. 



Kink Seeley Coefficients 




Kink Mass Shift Estimation 


n 


Cn{K) 


Nt 




1 


4.00000 








2 


2.66667 




2 


-0.663146 


3 


1.06667 




3 


-0.665798 


4 


0.304762 




4 


-0.666177 


5 


0.0677249 




5 


-0.666240 


6 


0.0123136 




6 


-0.666252 


7 


0.0018944 




7 


-0.666254 


8 


0.000252587 




8 


-0.666254 


9 


0.0000297161 




9 


-0.666255 


10 


3.12801 • 10"^ 




10 


-0.666255 




10 IS 



Table 1: Seeley coefficients (left), kink mass shift estimation (middle) and comparison between the shifts 
obtained from the modified (blue) and from the standard (red) asymptotic series (right) for several trun- 
cation orders in the (j)'^ model. 

In sum, the one-loop mass shift obtained from the modified heat trace expansion truncated at the 
order Nt = 10 is: 

A£;(<^K,A^i = 10) 



-0.666255 



13 



This estimation very precisely reproduces the exact result — ^ = — :^(3 — \/3arccos^) = — ^ = 
—0.666255 found in this model, see [8, 18]. 
Three points are worthy of mention: 

1. The Seeley coefficients for the 4>'^ kink obtained through the modified heat trace expansion are 
exactly the same as the Seeley coefficients for the sine-Gordon kink obtained from the standard 
heat trace asymptotics. 

2. It is clear from the above Figure that the convergence obtained upon using the modified approach 
is much better than in the standard procedure. The optimum asymptotic value is reached at a 
truncation order of three or four in the new method, whereas in the traditional approach more than 
ten terms were necessary. 

3. All this stresses how a better approximation to the K-heat function is obtained in the modified 
approach. The distortions introduced by the zero mode are almost suppressed when the whole 
effect of the kink fiuctuations is calibrated with respect to the simplest kink with a zero mode ! ! 



4.3 The cr = 3 and a = A parent potential models 

In reference [ a family of (l+l)-dimensional scalar field models governed by the action (1) with the 
following "parent" potentials^ 

orfi -I- 

u^'H<^) = ^l^{i-i-'m;h^]r (43) 

is fully analyzed. I~^{z; a,b) is the inverse incomplete regularized Beta function and the reason for the 
denomination "parent" is that the potential is identified by an inverse procedure. It is mandatory that 
the K-operator be 

-fv = — + a 



dx'^ cosh^x 
whereas the kink profile is identified from the zero mode of K: 

(^P{x) = sign(x) /[tanh^x; f] . (44) 

It is not completely clear from (43) and/or (44), but a = 1 and a = 2 respectively correspond to the 
sine-Gordon and X{4>)f_^_l models. Here we shall study the next two cases in depth (for a a positive 
integer): a = 3 and a = 4. 

In the first case, a = 3, the small vacuum and kink fiuctuation operators are 

Kq = — + 9 , K = — — r + 9-12 sech^x . 

Thus, we have v = 3 and V{x) = — 12sech^x and the normalized zero mode is: fo{x) = ^^^sech^x. 
The Seeley coefficients listed in the Table 2 are obtained from (34). We also specify the partial sums 
hK{P,Nt) (42) for distinct values of the truncation order Nt. The values of the one- loop mass shifts 
obtained from these partial sums are depicted in the attached figure by the blue line and compared with 
the data obtained in ' ]) using the standard procedure, which are represented by the red line. 
From the data in Table 2, the one-loop kink mass quantum correction is estimated to be 



AEi^ K,Nt = 10) 
^In it is explained that these models are well defined only for \cf)\ < 1 



-1.08450 

d 
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Kink Seeley Coefficients 




Kink Mass Shift Estimation 


Ti 




AT 




i 


1 9 nnnn 

iZ.UUUU 








z 


o/i nnnn 
Z4.UUUU 




z 


1 n^^i nQ 
— i.UoiUo 


KJ 








—1 07832 


4 


39.3143 




4 


-1.08262 


5 


34.7429 




5 


-1.08388 


6 


25.2306 




6 


-1.08429 


7 


16.5208 




7 


-1.08443 


8 


8.27770 




8 


-1.08448 


9 


3.89498 




9 


-1.08449 


10 


1.63998 




10 


-1.08450 




Standard Methfld 



Modified Method 



Table 2: Seeley coefficients (left), kink mass shift estimation (middle) and comparison between the shifts 
obtained from the modified (blue) and from the standard (red) asymptotic series (right) in the a = 3 
parent potential model. 



In [19] the exact kink mass quantum correction in this model is written in terms of the bound state 

— — — ^(6 — \/5arccos^ — 



^) = —1.08451. Therefore, the estimation obtained with the modified method is very precise. 



energies and the threshold of the continuous spectrum of the K operator: ^^^'^^^ 
\/8arccos 

In Figure 2a we plot the renormalized kink heat functions hxif^) for this model. In Figure 2b we depict 
the partial sums hK{P,Nt) (38) derived from the modified asymptotic series for several values of Nf. 
Finally, in Figure 2c we represent the partial sums hxif^jNt) (40), calculated by means of the GDW 
standard procedure. 



5 - 








1 2 





Figure 2: The exact function h^iP) (le-ft), partial sums /i_ft-(/3, Nt) from the modified factorization (middle) 
and partial sums /i_r-(/3, Nt) from the standard factorization (right) in the a = 3 parent potential model. 



In the second case a = 4, the small vacuum and kink fluctuation operators are 

d^ d^ 2 

Ko = — + 16 , K = — — J + 16 - 20 sech^x . 

Thus, we have u = 4 and V{x) = — 20sech^x and the normalized zero mode is: /o(x) = j^sech^x. 
In Table 3 we summarize the results obtained for this model showing the Seeley coefficients and the 
quantum mass shift estimations for different values of the truncation order Nt. 
The one-loop kink mass quantum correction is estimated to be 

AE{cl)K,Nt = U) 



-1.58003 . 



As in the previous example, we conclude that this result is very precise, recalling that the exact result 



is: 



:) 



^(10 



V3 



15 arcsin 



1 



V7 



arcsin ■ 



-1.58003, see [19]. 
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Kink Seeley Coefficients 




Kink Mass Shift Estimation 


Ti 


n ( i>r\ 

Cn\J^ ) 


AT 




i 


Z4.UUUU 








z 


yo.uuuu 




o 
z 


— i.oiiy / 


KJ 








—1 55773 


4 


705.829 




4 


-1.57144 


5 


1367.77 




5 


-1.57642 


6 


2206.55 




6 


-1.57843 


7 


3035.81 




7 


-1.57930 


8 


3632.62 




8 


-1.57969 


9 


3841.44 




9 


-1.57986 


10 


3637.21 




10 


-1.57995 



Modified Method • Stuidaid Method 



\ 



10 U 20 



Table 3: Seeley coefficients (left), kink mass shift estimation (middle) and comparison between the shifts 
obtained from the modified (blue) and from the standard (red) asymptotic series (right) in the a = 4 
parent potential model. 



Taken together, the results shown in the Tables 2 and 3 and represented in the above Figures confirm 
the pattern found in the X{4>)i+i model: (1) /ii<-(/5, Nt) matches the exact hxil^) at both small and large /3. 
hxiP, ^t), however, only behaves like the exact heat trace at high -g but, even though the approximation 

is better with increasing Nt, there is always a finitely large value of f3, such that hxi/S, Nt) starts to tend 
to instead of going to 1. (2) A much faster convergence to the optimum estimation is achieved using the 
modified GDW heat kernel factorization. (3) Some hierarchical structure arises: the Seeley coefficients 
for the kink heat trace modified expansion in the a parent potential model are those found in the a — 1 
model coming from the standard GDW factorization. This third point suggests that one could use any 
operator 

ax^ cosh X 

in the transparent Posch- Teller hierarchy as the starting point of the heat kernel expansion for any 
differential operator with a zero mode. 

4.4 A family of (</>)i+i models 

We now consider a one-parametric family of cp^ models. The potential, depending on the coupling constant 
a, is: 

[/(</.; a) = ^(</<2^a2)(c/>2- 1)2 , 
see Figure 3a. The vacuum orbit is Ai = {—1, 1} and the kinks that connect these vacuum points are: 



(t)K{x]a)= = . (45) 



In Figure 3b some members of this kink family have been depicted together with their energy densities. In 
these graphics it is clear that this model presents double kink solutions, a character accentuated when the 
parameter a approaches zero. Indeed, at the a — )• limit the potential reduces to U{(j), 0) = — 1)^, 

i.e., the pure (f)^ Lohe/Khare model [25]; a new vacuum point (/> = arises, and the vacuum orbit becomes: 
M. = {—1,0,1}. The \m\a^o(t)K{x]a) is a configuration where the (pxia) kink solution splits into two 
single kinks, one connecting the vacua (p = —1 and (p = and the other connecting = and (j) = 1. 
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Figure 3: Graphical representation of the potential (left), several kink profiles (middle) and kink energy densities 
(right) for a = 0.001, 0.01, 0.25, 0.5, 0.75, 1, 1.25, 1.5. 

The particle masses and the kink potential wells in this family of models are 

15(4a + l)2 6 (a2 + 3) (4a + 1) 



= 4(1 + 0^ 



V{x;a) 



2a cosh (2\/a2 + Ix) + 2a + 1 



2a cosh ( 2Va?~Tlx] + 2a + 1 



The Seeley coefficients of the modified asymptotic series evaluated using (34) have been implemented in 
the Mathematica code available in this paper for several selected values of a. Thus, the one-loop kink 
mass shifts are computed and are shown in Table 4 with a truncation of Nt terms. All these data are 
represented graphically in the attached figure in Table 4. 



a 


AE{cbK;Nt)/hjj 


Nt 


0.001 


-1.95321 


9 


0.01 


-1.65859 


9 


0.05 


-1.44956 


9 


0.1 


-1.35308 


9 


0.2 


-1.24166 


9 


0.3 


-1.17150 


10 


0.4 


-1.11564 


10 


0.5 


-1.07925 


11 


0.6 


-1.05903 


11 



a 


AE{cbK;Nt)/hjj 


Nt 


0.7 


-1.05493 


11 


0.8 


-1.06099 


11 


0.9 


-1.07728 


11 


1.0 


-1.10137 


11 


1.1 


-1.13185 


11 


1.2 


-1.16744 


11 


1.3 


-1.20727 


11 


1.4 


-1.25054 


10 


1.5 


-1.29667 


9 



Q.2 0.4 0.6 



1.0 1.2 1.4 ^ 



• ■ 



AEh 



Table 4: Quantum corrections to the kink masses in the cj)^ family for several values of the parameter 
a (left) and their graphical representation (right). The red squares represent the data obtained in the 
standard GDW approach, whereas the blue dots are the new data coming from the modified factorization. 

We observe in Table 4 that the quantum correction tends to infinity as a approaches zero. The reason 
for this stems from the fact that the vacua that are connected by the kink in the case a = give rise to 
different meson masses: 
this problem. 



[1] = 4 and §t3-[0] = 1, see [y] and references therein to see a treatment of 



4.5 The double sine-Gordon model family 

Finally, we shall apply the procedure to a family of models characterized by the action (1) where the 
potential is: 

C/((/);a) = 1 - (1 - a) cos(/) - acos(20) , 0<a<l,aGM . 

These models (here we are considering a as a free parameter) are referred to as the double sine-Gordon 
model in the literature, see [J-]. When a varies in the [0,1] range, U{(j);a) interpolates between the 
sine-Gordon, a = 0, and the re-scaled sine-Gordon, a = 1, potentials: 

Ui<P; 0) = 1 - cos = Ug""^ i^) , Ui<P; 1) = 1 - cos 2</> = U^^^ (0) . 
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If < a < 1, the vacuum orbit is Ai = {2Trk}, k £ Z. When a becomes greater than but is strictly less 
than 1, new relative minima (false vacua) emerge between each consecutive pair of vacua. At a = 1, the 
new minima become absolute minima, the vacuum points of the re-scaled sine-Gordon model, see Figure 
4 (left), such that the vacuum orbit becomes Ai = {irk}, /c G Z. The kink solution of each member of 
the family with a strictly less than 1 is: 



')K{x;a 



-2 arctan 



\/l — a sinh(-v/l + 3ax) 



TT + 2 arctan 



\/l — a sinh( vT+~3ax) 
VI + 3a 



(46) 



In Figure 4 we depict the kink profiles for several values of the coupling constant together with their 
energy densities. 






Figure 4: Graphical representation of the family of potentials for a — 0, 0.25, 0.5, 0.75, 1 (left), kink solution profiles 
(middle), kink energy densities (right) for a = 0,0.25,0.5,0.75,0.99,0.9999 in several double sine-Gordon models. 

For completeness we give the particle masses and the kink quantum potential wells: 

4(a - l)(3a + 1) ((15a + 1) cosh {2V3a + Ix) -9a + l) 



l + 3a 



V{x) 



-{a - 1) cosh {2^/3a + lx) +7a + l) 



We shall compute the one-loop quantum correction to the masses of the kinks (46) for several values of 
the parameter a within the modified asymptotic approach. These one-loop kink mass shifts are shown in 
Table 5, where the optimized truncation order is also depicted. Note that, depending on a, the asymptotic 
nature of the series fixes a lower or higher optimum number of terms to be considered in the formula 
(42). All these data are represented graphically in the attached Figure. 



a 


AE/hnid 


Nt 




a 


AE/hrud 


N 


0.0 


-0.318321 


20 




0.6 


-0.716205 


8 


0.1 


-0.370625 


3 




0.7 


-0.79149 


7 


0.2 


-0.429882 


5 




0.8 


-0.90928 


7 


0.3 


-0.495113 


5 




0.9 


-0.98757 


10 


0.4 


-0.565305 


6 




0.99 


-1.09687 


12 


0.5 


-0.637682 


7 




0.999 


-1.11267 


12 



Table 5: One-loop kink mass shifts in the family of double sine- Gordon models for several values of the 
parameter a (left) and their graphical representation (right). 

Bearing in mind that at the a — )• 1 limit the kink solution (46) is formed by two solitons of the re-scaled 
sine-Gordon model, it is tempting to interpret the quantum correction AE at this limit as the double of 



(n) 
sG 



the quantum correction of the sine-Gordon soliton. In particular, we recall that AE^ 
We observe behavior similar to this in Table 5, although the accumulation near the point a 
allow us to obtain the answer with complete precision. 



-0.636894;im. 
1 does not 
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5 Summary and outlook 



In a very brief summary and referring to outlooks we wish to stress again that we have proposed a 
modification of the Gilkey-DeWitt heat kernel expansion when addressing differential operators with 
non-trivial algebraic kernels. The new procedure has been implemented in several one-component scalar 
field theory models with both computational and conceptual gains with respect to the application in 
the same models of the standard procedure. We believe that the translation of these methods to more 
complicated models such as A^-component scalar field theory and/or Abelian Higss will be even more 
effective. 
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Appendix 

In this appendix we display a Mathematica code, which automatizes the computation of the quantum 
correction to the kink mass in one-component scalar field theory by applying formula (42), derived 
from the modified asymptotic series approach. The algorithm is divided into three subroutines: the 
identification of the density coefficients Cn{x,x) = ^Cn{x) by means of (34); the computation of the 
Seeley coefficients by integrating the density coefficients, finally evaluating the formula (42) to obtain an 
estimation of the kink mass quantum correction. 

• Calculation of the Cn{x,x) densities. 
The following Mathematica code 

densitycoef f icients [potential_, vacuuml_, vacuuin2_, lunaxj := Module [{varl , 
var2, var3, tomax, dl, v, vO, oper, fO, f6, x7, coef, k, coa, j, co} , 
(varl[phl_] = potential /. {y -> phl}; var2[phl_] = Simplify [PowerExpandCSqrt [2 
vaxl [phi] ] ] ] ; var3 [phl_] = Sign[var2 [(vacuuml + vacuuin2) /2] ] var2[phl]; coef 
= {}; v[x_] = Simplify [(D [varl [phi] , {phi, 2}]) /. {phi -> phi [x] }] ; vO 
= Sqrt [Simplify [(D [varl [phi] , {phi, 2}]) /. {phi -> vacuuml}]]; fO[x_] = 
(var3 [phi] / (Sqrt [Integrate [var3 [phi] , {phi, vacuuml, vacuum2}]])) /. {phi -> 
phl[x]}; dl[funj := Simplif y [(D [f un, x] ) /. {phi ' [x] -> var3 [phi [x] ] }] ; oper[fu8_, 
nl_] := Simplify [Nest [f 6, x7, nl] /. {f6 -> dl, x7 -> fu8}] ; tomax = 2 nmax; 
For[k = 0, k < tomax + 0.5, coa[0, k] = 0; k++] ; coa[0, 0] = 1; co[n_, kj : = 
Simplify [(l/(n + k) ) (coa[n - 1, k + 2] - Sum [Binomial [k, r5] oper[v[x] - v0"2, 
r5] coa[n - 1, k - r5], {r5, 0, k}] - 2 vO fO[x] oper[fO[x], k] KroneckerDelta[0, n 
- 1] - fO[x] oper[fO[x], k] (1 + 2 k) (2-n (v0)-(2 n - l))/((2 n - 1)!!))]; For[j = 
1, j < nmax + 0.5, tomax = tomax - 2; For [k = 0, k < tomax + 0.5, coa[j, k] = co [j , 
k] ; If[k == 0, coef = Append[coef, coa[j, 0]]]; k++] ; j++] ; Return[coef ] )] ; 

defines the module densitycoef f icients [potential., vacuuml., vacuum2_, nmax_] , which per- 
forms the work of calculating the coefficients c„(x, x). The arguments of this computational function 
are potential, the U{y) potential written by prescription as a function of the y variable, vacuuml 
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and vacuuin2, the two vacua connected by the kink solution in increasing order, and nmax, the 
A''-order truncation chosen in the computation of AE{(pK', ■ 

• Calculation of the Seeley coefficients. 

The Mathematica module seeleycoefficients [potential-, vacuuml_, vacuum2_, nmaxj de- 
pends on the same arguments as the previous one 

seeleycoef f icients [potential_, vacuuml_, vacuiim2_, nmaxJ := Module [{coef , densi, 
f, fl, f2, a = (coef = densitycoeff icients [potential, vacuuml, vacuum2, 

nmax]; fl[y_] = Simplify [PowerExpand [Sqrt [2 potential]]]; f2[y_] = Sign[fl [(vacuuml 
+ vacuum2)/2]] fl[y]; For[j = 1, j < nmax + 0.5, f [y_] = Simplify[(coef [[j]] /. 
{phl[x] -> y})/f2[y]]; a = Append[a, Integrate [f [y] , {y, vacuuml, vacuum2}]]; j++] ; 
Return [a] ) ] ; 

and provides us with the value of the Seeley coefficients. This subroutine calls the previous function 
in order to accomplishes its task. 

• Estimation of the quantum correction. 

The subroutine quantvmicorrection [potential., vacuuml., vacuum2_, nmaxJ 

quantumcorrection [potential-, vacuuml_, vacuum2_, nmaxJ := Module [{vO, corr, 
a}, (vO = Sqrt [Simplify [(D [potential, {y, 2}]) /. {y -> vacuuml}]]; a = 
Cliop[seeleycoeff icients [potential, vacuuml, vacuum2, nmax]]; corr = -(vO/Pi) - (1/(8. 
Pi)) (Sum[(a[[n]] (v0*(-2 n + 2)) Gamma[n -1]), {n, 2, nmax}]); Return [corr] )] ; 

completes the work by providing us with the one-loop kink mass quantum correction in the (1+1) 
dimensional scalar field theory model characterized by the potential term U{y). This function calls 
the two previous ones. 

The KinkMassQuantumCorrection_Modified.nb file containing this Mathematica code can be down- 
load at the web page http://campus.usal.es/~mpg/General/Mathematicatools, which includes examples 
and demos. We recommend this option in order to avoid transcription errors in the code. 
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